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Abstract 

In this work we investigate the issue of integrability in a classical model for non¬ 
interacting fermionic fields. This model is constructed via classical-quantum corre¬ 
spondence obtained from the semi-classical treatment of the quantum system. Our 
main finding is that the classical system, contrary to the quantum system, is not in- 
tegrable in general. Regarding this contrast it is clear that in general classical models 
for fermionic quantum systems have to be handled with care. Further numerical in¬ 
vestigation of the system showed that there may be islands of stability in the phase 
space. We also investigated a similar model that is used in theoretical chemistry and 
found this one to be most probably integrable, although also here the integrability is 
not assured by the quantum-classical correspondence principle. 
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1. Introduction 


The semi-classical approximation in quantum mechanics is by now an established 
tool to study in a comprehensible and pictorial, but still quantitative manner, non- 
classical characteristics of the quantum domain, in particular coherent phenomena due 
to quantum interference mm- 

The semi-classical program is in nature quite different from quasi-classical ap¬ 
proaches where the classical limit of the quantum theory is explicitly used in the form 
of trajectories that mark the evolution of expectation values of, for example, position 
and momentum. The regime of validity of the quasi-classical approximation is given 
by the possibility of reconstructing, at least in an approximate way, the quantum 
state using only its average position and momentum. A typical example is the time- 
dependent evolution of well localized, gaussian wavepackets [3]. Starting with the 
initial time, and assuming that the wavepacket essentially does not disperse, one can 
simply use the solutions of the initial value classical problem to guide the wavepacket’s 
evolution. It is apparent that the approximation fails when either the wavepacket can¬ 
not be adequately represented by a well localized Gaussian, or when different pieces of 
the initial distribution evolve in a non-rigid way and get to interfere with each other. 
In closed systems the time scale where these effects arise is called the Ehrenfest time, 
making reference to the short-time regime where the Ehrenfest theorem can be used 
to approximate quantum evolution by classical dynamics. 

The starting point of the semi-classical approach is, on the contrary, purely quantum- 
mechanical, namely, the path integral formulation [3]. One extremely important fea¬ 
ture of the path integral picture of quantum mechanics is that it allows us to define 
the classical limit of the theory. Contrary to the quasi-classical approach, the semi- 
classical methods do not require that the classical limit is in any way ’’correct”. Erom 
the point of view of the Feynman formulation, a classical limit is defined by the func¬ 
tional that appears as exponent in the integral over paths. 

The semi-classical approximation simply amounts for the evaluation of the path 
integral by the method of stationary phase (SPA) [T][^. It is of course satisfactory that 
the paths selected by the SPA are precisely the ones that satisfy Hamilton’s principle 
of least action, and therefore are ’’classical”. However, one must keep in mind that 
classical mechanics, from this point of view, is nothing but a set of rules that tell us 
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how to calculate actions that appear as exponents in semi-classical propagators. 


This tension between what one’s intuition indicates as classical limit and what 
the path integral demands to be used in order to construct semi-classical propagators 
reaches a maximum if the quantum theory we attempt to describe is a fermionic 
quantum held [B]. Here we must try to think classically about extremely counter¬ 
intuitive quantum manifestations of the kinematics of fermions, in particular the Pauli 
exclusion principle and the antisymmetry of the fermionic states. Nevertheless, a semi- 
classical approach to systems of interacting fermions is of utmost importance in several 
branches of physics, ranging from quantum chemistry to cosmology. 

Several approaches have been proposed to provide a classical theory that gives a 
glimpse of the dynamics of fermionic helds [7], but to our knowledge the only one 
that actually starts with the fermionic path integral and ends up with the associated 
van Vleck-Gutzwiller semi-classical propagator by strict application of the SPA was 
presented only recently [8] . The huge advantage of this approach is that one can be sure 
that, whatever strange or counter-intuitive features the classical limit has, it provides 
the unique consistent way to use a classical theory to study quantum interference 
effects for fermionic quantum helds through the systematic use of the semi-classical 
approximation. 

In this thesis, one of these apparently contradictory properties of the classical limit 
for fermionic helds is studied, namely the curious fact that, due to the intrinsic non¬ 
linearity of the saturation phenomena associated with the Pauli-principle, the classical 
limit of the free (non-interacting) fermionic held turns out to be non-integrable in 
general. This is in sharp contrast with the bosonic case [8], where non-interacting 
helds are paradigmatic examples of quantum systems with classical integrable limit. 

2. Outline of the quantum mechanical system 

Consider a system of L (bosonic or fermionic) spinless single particle states |i), i € 
{1 ,... ,L} (for example lattice sites) with eigenenergies ha and hoppings hij (where 
hij = hji to ensure hermiticity). Then the Hamiltonian of the system (in second 
quantization) is given by 


L 



( 2 . 1 ) 
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where a| and Oi are the creation and annihilation operators for the state i satisfying 



( 2 . 2 ) 


where [A,B]- := AB — BA is the commutator, which is used if the are of bosonic 
character and [A, B]+ := AB + BA is the anticommutator, used for fermions. 

As the matrix h is hermitian, we can find a unitary similarity matrix u to diag¬ 
onalize it. We then can use the matrix u to transform the creation and annihilation 
operators: 



(2.3) 



where we used the unitarity of u: (u’ = Uji. 

As the transformation is linear, the (anti-)commutation relations of the creation 
and annihilation operators are preserved: 


L 


L 



k,l=l 


k,l=l 


L 



(2.5) 


kj=l 


Using this transformation, H becomes diagonal: 


L 


L 


L 



( 2 . 6 ) 


where Wk = J2ij=i 'U'ikhijUjk are the eigenvalues of h. 


Because of eqs. (2.4l and (2.51, the Nk commute: 


[Nk,Ni] =AlAkAlAi-A]AiAiAk 

= A\.AiSki ± AlAjAkAi — {AjAkSki ± AjAlAiAk) = o (2.7) 


where “-I-” is used in case of bosons and ” in case of fermions. Because of this, the 
Nk obviously also commute with H = WkNk, thus we have found L independent 

constants of motion. This means the system is quantum-integrable. 
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In the original basis, the constants of motion read 

L 

Nk = A\Ak = ^ UjkUik&lcLj . ( 2 . 8 ) 

If the quantum-classical correspondence works well, we can expect the clasical system 
to behave in the same manners as the quantum one. Especially, it should also be 
integrable. 


3. Bosons 


The classical analogy in the bosonic case is very simple; the corresponding classical 
Hamiltonian for the bosonic case is obtained by replacing bi and b| by complex fields 
t/ji and ip* respectively [SlIS]. This may only be done for normally ordered products 
of creation and annihilation operators (which means all creation operators are to the 
left of all annihilation operators), so there is no ordering ambiguity. 

Given this, we can write down the classical Hamiltonian in the bosonic case: 

L 

(3-1) 

® j=i 


By the classical-quantum correspondence we expect conservation (of classical versions) 


of the quantities (2.8l: 


iVfe = ^ u^Uik'ipiipj 


(3.2) 


It has to be emphasized that these functions are just an educated guess for the con¬ 


stants of motion of the classical system (3.11. 

However, it can be shown by direct calculation, that this guess indeed gives L 
independent constants of motion, meaning that the classical system, in perfect analogy 
to the quantum system, is integrable. 


4. Fermions 

As in the bosonic case (section]^, the classical system is obtained by replacing the 
operators with complex functions. But because of the anticommutativity (versus the 
commutativity of the complex fields ip and ip*), the replacements are more complex. 
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In |8] they are given as 


c|ci(4.1) 

max(i,j) —1 

(l-2|V'fe|^). (4.2) 

fc=min(i,j) + l 

Here, however let us consider the slightly more general version 


c|ci->• , (4.3) 

max(i,j) —1 

-^i/)*V>j7(l^i|^)/(IV'jf) n (l-2|7fe|^) (4.4) 

A: = min(i,_ 7 ') + 1 

with an arbitrary smooth function f{x). 

Of course f{x) in reality is not completely arbitrary. For example setting f{x) = 0 
obviously would not be a reasonable choice. The limits of the freedom in the choice 
of / are investigated in m- The two most important choices are f{x) = e ^ and 
f{x) = a/1 - X. 

Due to these more complicated replacements, it is not possible (with affordable 
effort) to compute the Poisson-brackets for a Hamiltonian with a general number L 
of sites. So we start with the first non-trivial case, which is L = 2. 


4.1. Two sites 


With a fixed number of sites given, we can explicitly write down the matrix h deter¬ 
mining the Hamiltonian: 


h = 


A 


This form (with £ 1,62 € K) is required by the hermiticity of h. 


(4.5) 


Using the replacement rules (4.31 and (4.41 we get 


H =£i|V^i|" + £ 2 |V> 2 |" + (JV’l‘V '2 + J>l72‘)/(IV'1|')/(I72P) , (4.6) 

Nk = + \u 2 kf\lp 2 f + iuikU2k1pt‘lp2 + MlfeW 2 fe 7 lV’ 2 )/(IV'l|^)/(IV’ 2 |^) • (4.7) 


Given this it can be shown by explicit calculation, that the relevant Poisson-brackets 
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vanish: 


= Oyk G {1,2}, f G C\R) 


(4.8) 


and 

{Nk,Ni} = OWk,lG {1,2}, / G Ci(R). (4.9) 

This implies that for L — 2, the classical analogy to the fermionic quantum system is 
integrable. So now let us continue with L = 3. 


4.2. Three sites 


For the case L = 3 calculations similar to those in section 14.11 can be done. We 
will consider two different (yet very similar) systems: a linear one, where we only 
consider the interactions |1) |2) and |2) |3), and a circular system, where also 

the states |1) and |3) can interact. For simplicity we assume the same strength for all 
interactions. 

The linear Hamiltonian is given by the matrix 


h = 


ei J 0 

J* €2 J 

0 J* £3/ 


and for the cyclic sytem we get 


h = 


\ 


ei J J 
J* £2 J 

\r r £3y 


(4.10) 


(4.11) 


The calculations for the two cases are very similar and especially yield equivalent 


results, so here it shall suffice to discuss just the (shorter) one for eq. (4.101 in detail. 


4.2.1 Results 

The full calculations are very cumbersome, thus here we shall just state the result: 


{H,Nk} = 2i5l(^ - 2Jui*,U3fcV’iV'iV’3V’2/^(|V’i|^)/(IV’3|^)/(|V'2|^) 

- 2J*u{j,U3fcl/’lV'3^3V’2/(|t^l|^)/^(|V’3|^)/(|V’2|^) 
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+ JulkU3k1pl1p2 ((1 - 2|V>2|^)ff(|V’3|^) - 1 + SiV’al V^(IV’3|^))/(|V’l|^)/(l^2|^) 
+ J*UikU*3ki’3'lp2 ((1 - 2|V'2|^)ff(|V’l|^) - 1 + SiV’ll V^(IV’l|^))/(l^3|^)/(|^2|^) 
+ iJU3kU2ki’l1p3 + J*UikU2k'tpl1p3){9{\‘tp2f) “ (1 “ 2|V'2 | ^)) /(IV’l I ^)/(l V'3 |^)) 

(4.12) 


Now the question is: can we find a function f{x), such that this vanishes identically? 
The important part of the answer is: no. However, there are special cases (e.g. J = 0) 
where the Poisson-brackets do vanish, so we should investigate this issue in more 
detail. 

This detailed investigation is rather cumbersome and technical, thus here we shall 
just quote the result. The full considerations can be found in appendix [X| 

Ultimately, we find the following: the right hand side of eq. ( |4.12[ | can only vanish 
identically if 

JuifeM;,V’l>l>3V’2/'('I'?)/('k3)/('&2) = 0 . (4.13) 


However, reasonable assumptions (section A.ll state that every factor in this term is 
non-zero. 


4.2.2 Discussion 

In section[4.2.1|we stated that under reasonable assumptions (section[A.l[| the Poisson- 


brackets (4.121 cannot vanish identically. Thus there are two possibilities: either 
the Poisson-brackets really do not vanish or the assumptions were wrong. In this 
paragraph we will investigate this issue and see that both possibilities are partly true. 
In fact, the assumptions are legit for any generic setup, but there are special cases 
where we may not use them. In most of these special cases, it is very easy to see that 
the Poisson-brackets do vanish. 

The most important point is that the classical-quantum correspondence does not 
give full freedom of choice for /. Rather there are very restrictive conditions to be 
met, which essentially exclude pathological examples. Thus in sections ED and |A.2| 
we could innocently assume that everything is “nice”, e.g. there are points in phase 
space where the f’s and |V’|’s are non-zero. 

To be precise, the most important example is f{x) = e~^ which certainly fulfills 
the non-zero conditions by vanishing nowhere in the phase space. Thus the proof 
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certainly holds for this case. 

Especially in the field of theoretical chemistry, f(x) = \/l — x is widely used. This 
provides certain particularities. First of all, the domain of / is bounded to a; < 1. 
Thus we now have a boundary of the domain, which we have to consider separately. 
Assume we start with the system in a configuration where = 1 Vi. As shown below 


(eqs. (4.141 and (4.151), the quantity IV’iP is ^ conserved quantity. But because of 
the condition I'l/’iP < 1, this means that the I'i/’iP’s can neither increase nor decrease, 
because as one of them got smaller, another one would have to grow, which is not 
possible. Thus ~ const. = IVi, which directly implies /(IV'iP) = const. = OVi. 
Thus the Poisson-brackets do vanish on any trajectory that runs on the border of the 
domain. 

However, besides this special case all the conditions also apply for f{x) = — x, 

so inside the domain the Poisson-brackets still do not vanish. 


4.2.3 Consequences 


Now we shall discuss the implications of the fact that the Poisson-brackets (4.121 


do not vanish. The bottom line is that for L = 3, the Nk are in general no longer 
constants of motion. This does not directly imply that the system is not integrable 
anymore, but our educated guess for the constants of motion is just wrong. 

But still the Nk are not completely useless. The interpretation of Nk in the quan¬ 
tum system is the number of particles in state k. Of course, the overall number of 
particles N = in the system is conserved, and this stays true in the classical 

analogon. To show this, we note that all the terms in eq. ( |4.12[ ) are proportional to 
a combination of Uij’s of the form UikU*k where i yf j. Employing the unitarity of u, 
these expressions, when summed over k, yield Sij which is zero for i yf j. Thus 

3 

{H,iV} = ^{H-,iV4 = 0, (4.14) 

k=l 

SO the total particle number, as stated above, still is conserved. It can be shown that 


N = \^l\^ -I- \lp2\^ + \lp3\^ ■ 


(4.15) 


The last thing to mention about eq. (|4.12[| is that all terms are proportional to J 
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or J*, so a system of three independent fermionic states (which means J = 0) still is 
integrable, which is also trivially seen from its separability. 


5. Numerical investigation of the fermionic L=3 system 

From the results of section fF^ we can make no statement concerning the integrability 
of the system in general. To investigate this issue, we simulated the system and looked 
at Poincare-surfaces of section. In order to do so, we first transform the canonical 
variables in a suitable way which is discussed in section [57l| The numerical results are 
shown in section IFTil 


5.1. Variable transformations 


In the above considerations we saw that for L = 3 the direct analogy between the 
quantum and the classical system does not give enough constants of motion to ensure 
integrability anymore. However, there still are constants of motion: of course H itself. 


as it is time-independent, is a constant of motion and in eqs. (4.141 and (4.151 we saw 
that the total particle number 


N = Ni + N2+N3 = IV'il" + |V'2|" + IV'al" 


(5.1) 


also is a constant of motion. Using this, we can reduce the dimensionality of our 
problem by variable transformations. 

First of all, we want to get rid of the complex variables tpi and and instead 
express our system by means of real variables. Doing so we have to take care about 
the fact that the equations of motion for i/ji and 'tpt (which are given for example 
in [8] eqs. (4.2a,b)) are not the standard Hamilton’s equations. Instead there is an 
additional factor i in them (which is needed to assure consistency). But for the real 
variables we are about to introduce we must have Hamilton’s equations without this 
prefactor. Thus the first transformation introduced here cannot be canonical, as it 
will not precisely preserve the equations of motion. 

However, we can get close to the needed transformation by considering the gener- 
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ating function F = X]i=i fV'?® which yields the (canonical) transformation 

,* dF _ 2 i 0 _ dF . 2 - 2 ie ■ i ,* /c r,\ 

Vi = = V’e , ^ ® = -WiVi ■ (5.2) 

Now we introduce the coordinates rii ~ihi = ■ This cancels the i in the equations 

of motion, so in the real variables we obtain Hamilton’s equations 


rii 


dH ■ dH 

de ' dm' 


(5.3) 


From now on we will consider the h given in eq. (4.11 1 . Using the new coordinates rii 
and 6 i, the Hamiltonian then reads 


FI = mei + 712^2 + Uses 

+ (Je-‘®i+‘®" + J*e‘®i-‘®")yriT7i^/(ni)/(n2) 

+ (Je“‘®^+‘®’^ + J*e‘®’^-‘®’^)yn^/(n3)/(n2) 

+ (Je“‘®^+‘®’^ + J*e‘®^“‘®’^)(l - 2n2)Vnin3 /(ni)/(n3). (5.4) 


Note that eq. (4.151 in these variables reads 


N = ni + 712 + n 3 . (5.5) 

Now we will exploit the conservation of N to reduce the dimensionality of the problem. 
Therefore we again introduce new coordinates ((n, a), (m, j3), [N, 0)) where we want 
N to be the total particle number. We transform to these coordinates by means of 
the generating function 


F{n,m, N, 0i, 6 * 2 , 6 * 3 ) = ?r(^'i — 82 ) + 771(63 — 62 ) + N 92 , (5.6) 


which yields the transformations 

<^ = 1^=61-62, ^ = ^=63-62, e = 

dF dF \T dF 

ni = f^=n, n2 = f^=Ar-ni-n3, n3 = |^=m. 
This trivially gives N = ni + 712 + 113 which is what we intended to do. 


(5.7) 
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Rewriting the Hamiltonian (5.41 in these coordinates, we obtain 


H = n(ei — £2) + Ne2 + m{t3 — £2) 

+ (Je~‘“ + J*e^°‘)y/n\/N — n — m f{n)f{N — n — m) 

+ — n — m f{m)f{N — n — m) 

+ (Je”"^”*"'^ + J*e'“~‘^)(l — 2{N — n — m))\/nm f(n)f{m) . (5.8) 


Note that H does not depend on 0. So again we see that N, as the conjugated 
momentum to the cyclic coordinate 0, is conserved. In fact this Hamiltonian can be 
interpreted as a system with two degrees of freedom (instead of the previous three), 
where N is just a constant parameter. 

The variable pairs (n, a) and (m, /3) still are of the nature of polar coordinates. 
For visualization (see below) we want to use coordinates of cartesian nature, thus we 
perform one last transformation. Unfortunately, the transformation to coordinates 
{x 3 ,X 2 ,yi, t/ 2 ) where x\ X 2 = n and yl + y 2 = m ]s not canonical. Nevertheless we 
will use these coordinates, which is justihed by the following argumentation: 

Consider the generating function 

I o 1 o 

F = -ail tana +-i)i tan/3 . (5.9) 

This ultimately gives the transformations 


^ cos a , ^ sin a 

^ = y/m cos /3 , ^ = sin P . 


(5.10) 


so obviously Vi = However, if we use these coordinates, the equations of 

motion are no longer the canonical ones. Consider for example the canonical equation 
of motion for ii (for brevity we write {qpi = 1,2, 3, 4} = {xi,X 2 ,y\,y 2 })’- 


dH _yX dH dq^ _y^dH 5i2 _ 1 dH 
dx2 ^ dqi dx2 ^ dqi y/2 y/2 dx2 ’ 

1=1 Z=1 * * 


(5.11) 


With xi — this gives Xi = Analogous we obtain the factor | in every 

equation of motion. However, this is the only place where the non-canonicity of the 
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coordinate transformation comes into play, and on the other hand there is no other 
place in the whole system, where the time coordinate appears. Thus we can absorb the 
constant factor into onr definition of time and thns re-obtain the canonical equations 
of motion. In the end, this just means that the evolution of the system is twice as 
slow as it would be otherwise, but (as the Hamiltonian is time-independent) the overall 
speed of the evolution of the system is irrelevant for the evolution itself, which is what 
we are interested in. 

So finally we write the Hamiltonian as 

H = {xi + X 2 ){ei — £ 2 ) + ( 2/1 + y2)i^3 — £ 2 ) + Ne 2 

+ 2^{J{xi - 1 x 2 ))^N -xl-xl-yl-yl f{x\ + xl)f{N -xl-xl-yl- yl) 
+ 2K(J(2/i + m))\jN -xl-xl-yl-yl f{xl +xl)f{N -xl - xl - yl - yl) 

+ 2K(J(xi - ix 2 ){yi -I- ij/ 2 ))(l - 2(iV - xl - xl - yl - yl))f{xl + xl)f{yl + yl). 

(5.12) 


5.2. Poincare-surfaces of section the fermionic three-site system 

In this section we require some general knowledge of theoretical mechanics and the 
technique of Poincare-surfaces of section in general. An introduction to this topic can 
be found for example in in. 

In order to do a Poincare-surface of section, an important step is to solve the 
eqnation H = E for one of the canonical variables. Unfortunately, in onr specific con¬ 


siderations, the Hamiltonian (5.8 1 does not have a form which admits a nice solution 
of this equation for any of the variables. Thus we cannot easily find a parametrization 
of the energy shell. But we can use an alternate procedure: first we slice the phase 
space by looking at the coordinate hyperplanes of a canonical variable (now we see 
why we wanted to transform to coordinates of cartesian character: this better suits the 
intuitive picture of R") and then we manually restrict the image to the energy shell. 
This restriction is performed by considering only trajectories with a given energy. 

To get familiar with the energy shell of our concrete problem, an example is de¬ 
picted in fig. HD 
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Figure 5.1: Several slices of the energy shell of the Hamiltonian (5.8) for a generic choice 
of parameters. These images are inteded to give an overview over the full (4-dimensional) 
energy shell. Therefore we fixed the canonical variable j/i, thus the three dimensions of 
every small image represent the remaining three canonical variables. The values for yi 
increase from top-left to bottom-right in non-equidistant steps. 


5.3. Interpreting Poincare-surfaces of section 

In this work, the main reason for doing Poincare-surfaces of section is to distinguish 
between integrable and non-integrable systems. This is very easy, given the following 
considerations: the system we are looking at has two degrees of freedom and the 
Hamiltonian is a constant of motion. Thus the question of integrability reduces to 
the question whether there is another constant of motion or not. If there was a 
second conserved quantity, the invariant manifolds of the flow of the system would be 
two-dimensional, thus each of them would appear one-dimensional in the Poincare- 
surface of section. If, however, there was no second constant of motion, each invariant 
manifold would be three-dimensional and thus its image in the Poincare-surface of 
section would be two-dimensional. 

Of course the numerics do not directly give pictures of the invariant manifolds of 
the system, but just distinct points of single trajectories, where each trajectory runs 
on an invariant manifold. But if the trajectory covers a significant part of the invariant 
manifold (especially if the trajectory is a dense subset of the manifold), we get a good 
image of what this manifold looks like. Of course the picture gets better, the longer 
the trajectories are traced. 

So finally the prescription for interpreting the shown Poincare-surfaces of section 
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is as follows: if the points of one trajectory form a (one-dimensional) line, there is a 
second constant of motion and the system is integrable, whereas if the points spread 
over a certain (two-dimensional) area, the system is not integrable. 


5.4. Numerical results 

For the numerical calculations we considered the cyclic system (see section [F^ and eq. ( |4.11[ |). 
The computations where performed using Wolfram’s Mathematica 10. 

The main results of the numerical calculations can be summarized in two pictures. 


which are fig. 5.2 for the case f{x) = e ^ and fig. 5.3 for the case f{x) = — x. The 

implications of these results are explained in the figure captions respectively. 
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Figure 5.2: Poincare sections of several trajectories (marked by color) of the system (5.12) 
with f{x) = at Cl = £2 = £3 = 1, J = 0.6, N = 3, H = 3.14. Apparently there are 
regions on the energy shell (which is depicted in yellow), where the system displays non- 
integrable behaviour, whereas it looks (close to) integrable in other regions. The presence of 
the non-integrable regions proves the result of the calculations in section 4.2 where we saw 
that the educated guess for global constants of motion does not work. Here we see that this 
is because the system simply is not integrable globally. This implies a partial breakdown 
of the classical-quantum correspondence. However the “close to integrable” regions imply 
the existence of a local second constant of motion. This is under current investigation. 
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Figure 5.3: Generic Poincare-surface of section of the system (5.12) with f{x) = y/1 — x. 
According to section ^ this definitely looks like there should be a second constant of 
motion in this case. However, direct calculation (section 4.21 shows that the classical- 
quantum correspondence investigated in this work still does not give enough constants of 
motion. Thus we conjecture that there is an additional constant of motion which we have 
no knowledge of. 
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6. Conclusion 


At the beginning of this work it was shown by direct calculation, that the quantum 
system under consideration is integrable both in the bosonic and the fermionic case 
and for arbitrary number of sites. 

Using the quantum-classical correspondence from [S], we thus expected the classical 
Hamiltonian to be integrable. Direct calculations proved that this is the case for the 
bosonic system, for the fermionic system, however, it was shown that the classical- 
quantum correspondence in general does not give enough constants of motion. 

To further investigate the system, it was simulated numerically, which resulted in 
the insight that the fermionic three-site system is not integrable in general. This breaks 
the classical-quantum correspondence. The reasons for this behaviour are subject to 
current research. 

However, the numerics also showed integrable regions in the Poincare-surfaces of 
section, which indicates an additional (local) constant of motion which we have no 
knowledge of. In the case of f{x) = Vl — x there even seems to be another global 
constant of motion which would imply that this system is integrable. However, also 
for f{x) = we could show that the quantities we guessed from the classical- 

quantum correspondence are no longer conserved. 
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Appendices 


A. Investigation of expression (4.12) 


A.l. Technicalities 


As they only depend on derivatives of the analytic functions H and Nk, the Poisson- 
brackets are analytic functions in the canonical variables tpi and ip*. For an analytic 
function, to vanish identically is equivalent to the function and all its (partial) deriva¬ 
tives vanishing at some point. So in order to prove that the Poisson-brackets do not 
vanish identically, it is sufficient to find a point where an arbitrary partial derivative 
does not vanish. 

To exclude trivial solutions, assume that JuikUjk yf 0 for arbitrary i yf j. Now 
take a point in phase space where /(|V’i|) 7 ^ 0''^* and \ipi\ yf OVz. It is reasonable to 
assume f{x) yf 0 on the inside of the domain of /, thus such points do exist on any 
trajectory that does not solely run on the border of the domain of /. 

Now we have to find suitable partial derivatives. Therefore note that if we define 
Ip = 4'e'* (4^, $ G R), we get 


dip dip* 

df{ip,ip*) ^ df{ip,ip*) _ df{ip,ip*) 
dip dip* 

Thus we can write the Wirtinger derivatives as 


d e-'* 

( ^ 

i d\ 

d e'* 

( ^ + 

i d \ 

dip~ 2 ' 


vE' 5$/ ’ 

dip* ~ 2 


in d^J 


From this we see that 


df{ip,ip*) ^ ^ df{ip,ip*) df{ip,ip*) ^ df{ip,ip*) 

dip dip* d^ d^ 


(A.l) 

(A.2) 

(A.3) 

(A.4) 


Thus it will suffice to look at derivatives with respect to the phases of the ipi. This 
drastically simplifies the problem, because / and g only depend on |^i| = 4^^. 

Now we need to know how to handle the ^-symbol in the formula. By definition 
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^/= ^(/-/*). so 


^ d^f df{ip,ip*) d^f df*{ip,ip*) 

9$ df 9$ 9/* 9<& 

^ 1 9/(-i/>,V)*) _ 1 drill), j}*) ^ 9/(V),V)*) 

2i 9$ 2i 9$ 9$ 

Here we used that 

9/*(v>,v<*) ^ /1 9/(v>,v<*) y ^ /1 9/(V), v>*) y 
9$ V2i 9$* y V2i 9'1> / 

because $ £ R. Finally we can state that ^ and commute. 

For briefness in the calculation also note that 


9’"i/>r 

W]r 




9™(V-*)^ 

(9$,)™ 




and 9'(i/) = IR/. 


(A.5) 


(A.6) 


(A.7) 


A.2. The derivatives of the Poisson-brackets 

Now we will look at different derivatives of {H, N^} and impose they all vanish. This 
will lead to a contradiction with our assumptions. 

In the computations we will implicitly use the relations of section |A.1[ but as it 
would unnecessarily elongate this section, we will not explicitly point out every single 
step. 

+ (A.8) 

= -KJ*utj,U3fcl/>lV)3y’3V’2/(5'?)/^(4'3)/(VE'2) ■ 

(A.9) 

+ J*wtfcW3/cy’lV>3>3>2/('I'?)/('I/l)/('k2)) (A.IO) 


21 















aJuifcn;,V'l>l>3V>2/"('&?)/('I'3)/('l>2) = -aJ*uI,M3feV'l^3>3>2/(«'?)/"('I>l)/('I'^). 

(A.ll) 

° = = 8ia( - 2 Jui,«;,V’iV'i>3V>2/"('&?)/('I>3)/('I'2) 

(a$i) a$39$2 ^ 

+ J*Mtfe«3.^lV>3>3>2/('I'?)/"('I'i)/('I'^)) (A.12) 


29JuifeM3fc?/’l*V'r^3V’2/^(^'?)/(5'3)/(^2) = ‘^J*'ulkU3k^li’3'tp3^2fi'i’1)fi'i’l)f{'i'l) ■ 

(A.13) 


0 = '2 = 8i5R( - 2Ju,ku;krM3^^2f{■^^)f{<^>^)f{<^>^) 

(9$l) 9$3 (9$2) ^ 

+ J*utfeU3fc'!/’lV’3V’3^2/(^'?)/^(^'3)/(5'2)) (A.14) 

2KJui,«LV’t^l>3V'2/"(®?)/('I'3)/('I'2) = KJ‘^itfcn3fe^lV'3>3>2/('I'?)/"(®3)/('I'2) • 

(A.15) 


From eqs. (A.9l, (A.ll I, (A.131 and (A.151 we can conclude that 


KJuifeu5fcl/jJ'V'l^3V>2/^(5'l)/(4/3)/(’I'i) = '^Juikulk'>pl‘4>*l'ill3'>p2f{^\)f{^3)f(^2) = 0 

(A.16) 

and thus 

JuifeU^fe^)>l>t^3V’2/"(4'?)/(^3)/(^i) = 0 . (A.17) 

This however, is in direct contradiction to our assumptions, because basically we 
assumed every part of this expression to be non-zero in section[A.l| 
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